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Abstract 

The method of positive commutators, developed for zero temperature prob- 
lems over the last twenty years, has been an essential tool in the spectral analy- 
sis of Hamiltonians in quantum mechanics. We extend this method to positive 
temperatures, i.e. to non-equilibrium quantum statistical mechanics. 

We use the positive commutator technique to give an alternative proof of 
a fundamental property of a certain class of large quantum systems, called 
Return to Equilibrium. This property says that equilibrium states are (asymp- 
totically) stable: if a system is slightly perturbed from its equilibrium state, 
then it converges back to that equilibrium state as time goes to infinity. 

Keywords: positive commutator, Mourre estimate, return to equilibrium, virial theorem, 
Fermi golden rule 
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1 Introduction 



In this paper, we study a class of open quantum systems consisting of two inter- 
acting subsystems: a finite system, called the particle system coupled to a reservoir 
(heat bath), described by the spatially infinitely extended photon-field (a massless 
Bose field). The dynamics of the coupled system on the von Neumann algebra of 
observables is generated by a Liouville operator, also called Liouvillian or thermal 
Hamiltonian, acting on a positive temperature Hilbert space. Many key properties 
of the system, such as return to equilibrium (RTE), i.e. asymptotic stability of the 
equilibrium state, can be expressed in terms of the spectral characteristics of this 
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operator. 

Applying the positive commutator (PC) method to the Liouville operator of sys- 
tems in question, we obtain rather detailed information on the spectrum of these 
operators. This allows us to recover, with a partial improvement, a recent fundamen- 
tal result by several authors on RTE. 

Our main technical result is a positive commutator estimate (also called a Mourre 
estimate) for the Liouville operator. This result holds for a wider class of systems 
than previously considered. 

Spectral information on the Liouville operator, and hence the property of RTE, is 
extracted from the PC estimate through Virial Theorem type arguments. It turns out 
that the existing Virial Theorem techniques are too restrictive to apply to positive 
temperature systems, and we need to extend them beyond their traditional range of 
application. 

There is a restriction on the class of systems for which we prove RTE, due to our 
Virial Theorem type result mentioned above. This is the first result of this kind, and 
we expect that it will be improved to yield the RTE result for a considerably wider 
class of systems. 

1.1 A class of open quantum systems 

The choice of the class of systems we analyze is motivated by the quantum mechanical 
models of nonrelativistic matter coupled to the radiation field, or matter interacting 
with a phonon field (quantized modes of a lattice), or a generalized spin-boson sys- 
tem. For notational convenience, we consider only scalar Bosons. A good review of 
physical models leading to the class of Hamiltonians considered here is found in [HSp] . 

The non-interacting system. The algebra of observables of the uncoupled 
system is the C*-algebra 21 = B(H P ) <g> W($)o), where B{H P ) denotes the bounded 
operators on the particle Hilbert space H p and W(f)o) is the Weyl CCR algebra over 
the one-particle space $) = {f <E L 2 (R 3 ,d 3 k) : j \k\~ l \f{k)\ 2 < oo}. The restriction 
to / G S)o comes from the fact that we will work in the Araki- Woods representation 
of the CCR algebra, which is only defined for Weyl operators W(f) with / e f)o 
(see [AW], [JP1,2], [BFS4]). The dynamics of the non-interacting system is given 
by the automorphism group K. 3 t i— > a ti o € Aut(2l), a t fl(A) = e ttH ° Ae~ ltH ° , where 
H = Hp <g) If + lp <g) Hf is the sum of the particle and free field Hamiltonians. 

Hq acts on the Hilbert space H p <g) Hf, where Hf = 0^L o ^o sym * s the Fock space 
over Sjo and Hf is the free field Hamiltonian, i.e. the second quantization of the 
multiplication operator by u = \k\, Hf = dr(u;); if a*(k), a(k) denote the (distribution 
valued) creation and annihilation operators, then we can express it equivalently as 
Hf = J uj(k)a*(k)a(k)d 3 k. The particle Hamiltonian is assumed to be a selfadjoint 
operator on H p which has purely discrete spectrum: 

a(H p ) = {£,}°° =0 , (1) 
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(where multiplicities are included, i.e. for a degenerate eigenvalue Ei, we have Ei = Ej 
for some j ^ i), and we denote the orthonormal basis diagonalizing H p by {ifij}. Let 
tr denote the trace on B(7i p ), then we further assume that 

Z p {(3) := tie-P** < oo, V/3 > 0. (2) 

We do not need to further specify the particle system. As a concrete example, one 
may think of a system of finitely many Schrddinger particles in a box (hence the 
name particle system), or a spin system. In some of our results (see Theorem 4.4 on 
the Fermi Golden Rule Condition), we shall assume that the spectrum of H p is finite 
(N- level system). 

The equilibrium state at temperature T = 1/(3 > for the non-interacting system 

is given by the product ojp$ — uj P p®ujp G 21*. Here, Wp(-) = ^rp^g; - is the particle- 

Gibbs state at temperature (3 and Up is the field /3-KMS state that describes the 
infinitely extended field in the state of black body radiation, i.e. its two-point function 
is given according to Planck's law by ajt(a*(k)a(k')) = ) _ ^he QNg construction 
for (A, cx t fl, oJp,o) yields the (up to unitary equivalence) unique data (TC, L , flpfi, ir) 
(dependent on (3). Here, TL is the GNS Hilbert space with inner product (•, •), Qp j0 
is a cyclic vector for the *-morphism it : 21 — > B(H) (the representation map), and 
the Liouvillian Lq is the selfadjoint operator on 7i implementing the dynamics, i.e. 
satisfying L Qp^ = and 

up, (a t ,o(A)) = (Qp, , e itL °7T(A)e- itLo Qp,o) , VA e 21. 

This GNS construction has been carried out in [AW] (for the field, the particle part is 
standard since it is a finite system), see also [JP1,2], [BFS4]. We shall not explicitly 
use the representation map 7r here and thus omit its presentation which can be found 
in the above references. The GNS Hilbert space and cyclic vector are given by 

H = H p ®Ti p ®F{L 2 (^LxS 2 )), (3) 

fy >0 = (4) 

where is the particle Gibbs state at temperature (3 given in (21). JF(L 2 (IR x S 2 )) 
is the Fock space over L 2 (R x S 2 ) with vacuum f2, which we call the Jaksic-Pillet 
glued space. It was introduced by Jaksic and Pillet in [JP1] and is isomorphic to 
Hf <8> Hf, the field GNS Hilbert space constructed in [AW]. It is easily verified that 
the Liouvillian is given by L = L p + Lf (see also [JP1,2]). We write simply L p instead 
of L p (g) $-f(l 2 (RxS 2 )) an d similarly for Lf. Here, L p = H p ® l p — l p ® H p , Lf = dr(u) 
and u is the first (the radial) variable in 1 x S 2 . It is clear that the spectrum of L p 
is the discrete set {e = Ei — Ej : E it j G cr(H p )} and the spectrum of Lf is the entire 
real axis (continuous spectrum) with an embedded eigenvalue at (corresponding to 
the vacuum eigenvector Q). Consequently, Lq has continuous spectrum covering the 
whole real line and embedded eigenvalues given by the eigenvalues of L p . 
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The interacting system. We now describe the interacting system by defining 
an interacting Hamiltonian acting on H p <S> H-f 



H = H + Xv, 



(5) 



where the coupling constant A is a small real number, and 



v = G®(a(g) + a*(g)). 



(6) 



Here, G is a bounded selfadjoint operator on H p . The function g e 5)o is called the 
form factor and the smoothed out creator is given by a*(g) = J d 3 k g(k)a*(k). We 
assume g to be a bounded C 1 -function, satisfying the following infra-red (IR) and 
ultra-violet (UV) conditions (recall that uo = \k\): 



In addition, we assume that conditions (7) hold for the derivative d^g, if p, q are 
replaced by p — 1, q + 1. 

We point out that the value coming from the model of an atom coupled to the 
radiation field in the dipole approximation is p = 1/2 (without this approximation, 
p = —1/2). From now on we will refer to p = 1/2 as the physical case. 

The interacting Hamiltonian (which describes the coupled system at zero temper- 
ature) corresponds to an interacting Liouvillian (positive temperature Hamiltonian) 
which is given by (c.f. [JP1,2], [BFS4]): 



Here, Gi := G (g> i p , G r := l p ® CGC, where C is the antilinear map on H p that, in the 
basis that diagonalizes H p , has the effect of complex conjugation of coordinates. The 
origin of C is the identification of the Hilbert-Schmidt operators on 7i p with 7i p ® 7i p 
via the isomorphism <-> y 9 ® Cip (see also [JP2], [BFS4]). Moreover, we have 



and g2(u,a) = —gi(—u,a), where the function \x = fi(k) is the momentum density 
distribution, given by Planck's law describing black body radiation: ji(k) = (e^ w — 
1) _1 , uo = \k\. The structure of g\ in (10) comes from the Jaksic-Pillet gluing which 
identifies L 2 (IR 3 ) ©L 2 (IR 3 ) with L 2 (R + x S* 1 ) via the isometric isomorphism (f\, f 2 ) h- > 
/, f(u, a) = ufi(u, a) for u > and f(u, a) = uf 2 (—u, a) for u < 0. For more detail, 



IR: \g{k)\ < Cu p , for some p > 0, as uj — > 0, 

for some results, we assume p > 2, 
UV: \g{k)\ < Coj~ q , for some q > 5/2, as uj — > oo. 



(7) 



L = L + XI, 

I = G l ®(a*(g 1 ) + a(g 1 ))-G r ®(a*(g 2 ) + a(g 2 )). 



(8) 
(9) 



defined, for g G L 2 (M + x 5 2 




(10) 
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we refer to [JP1,2]. 

For A 7^ 0, one can construct a vector ftp ^ G H s.t. the state defined by ujp^A) = 
(flp,\, AQp^) is a (3-KMS state w.r.t. the coupled dynamics a t (A) = e ttL Ae~ ltL , where 
A is an element in the von Neumann algebra 971 := B(H P ) <S> 7r(W($)o)) ( we& k closure 
in i3(jF(L 2 (M x S 2 ))) ). An extension of the algebra of observables to this weak closure 
is necessary since the full dynamics does not leave B(J~t p ) <8> 7r(W(^o)) invariant. It is 
not difficult to show that (9Jt, a t ) is a W*-dynamical system (compare also to [FNV], 
[JP2]). Notice in particular that LQp X = 0. 

The construction of Qp : \ goes under the name structural stability of KMS states, 
see [BFS4] for this specific model, but also [A], [FNV], [BRII]. For (5\X\ small, one 
has the estimate (for the 0-notation, see after (20)): 



We show in Appendix A.l that L is essentially selfadjoint (Theorem A. 2). 
1.2 Spectral characterization of RTE 

We define the equilibrium states at temperature T = l/,3>0tobe the /3-KMS 
states. Hence the equilibrium state of the coupled system at inverse temperature 
(3 > is given by the above constructed ujp t \ G 9JT*. A conjectured property of KMS 
states is their dynamical stability (which should be a natural property of equilibrium 
states). In our case, this means that uj' o a t — ■> up t \ as t — > oo, for states uj' that 
are close to u>p,\- This is called the property of return to equilibrium. Apart from 
specifying the mode of convergence, it remains to say what we mean by uj' is close to 
ojp,\- There is a natural neighbourhood of states around uip y \ in which the dynamics 
is also determined by L: the set of all normal states uj' w.r.t. 0Jp,\- By definition, uj' 
is normal w.r.t. ojp,\, iff 



where tr(-) is the trace on the GNS Hilbert space H given in (3), p is a trace class 
operator on Tt, normalized as trp = 1. 

Proposition 1.1 (spectral characterization of RTE). Let 3Jt C B(H) be a 

von Neumann algebra and suppose that ujp(-) = (dp, ■ ftp) : 9JT — > C is a f3-KMS state 
with respect to the dynamics a t G Aut(djX). Suppose that the Liouvillian L generating 
the dynamics on H has no eigenvalues except for a simple one at zero, so that the 
only eigenvector of L is Vtp. Then, for any normal state uj' w.r.t. ujp, and for any 
observable A G 9JI, we have 



np,x-np, \\ = o(p\\\)- 



(ii) 



VAeM: uj'(A) — tr(pA), 



(12) 




(13) 
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This means that the system exhibits return to equilibrium in an ergodic mean sense. 

The proof is given e.g. in [JP2], [BFS4], [M]. Better information on the spectrum 
of L yields stronger convergence; if L has absolutely continuous spectrum, except a 
simple eigenvalue at 0, then (13) can be replaced by lim^oo u'(a t (A)) = <jjp(A). 



1.3 The PC method 

This section introduces the general idea of the PC method. As we have seen above, 
the Liouville operators in the class of systems we consider consist of two parts: 

L = L + XI, 

where L is the uncoupled Liouville operator, describing the two subsystems (particles 
and field) when they do not interact. / is the interaction, and A is a real (small) cou- 
pling parameter. The spectrum of L consists of a continuum covering the whole real 
axis, and it has embedded eigenvalues, arranged symmetrically w.r.t. zero. Moreover, 
zero is a degenerate eigenvalue. We would like to show that for A ^ 0, the spectrum 
of L has no eigenvalues, except for a simple one at zero, because then Proposition 1.1 
tells us that the system exhibits RTE! 

x ^ m 



XX X XX 



o o 

degenerate non-degenerate 



In other words, we want to show that all nonzero eigenvalues of L are unstable 
under the perturbation XI, and that this perturbation removes the degeneracy of 
the zero eigenvalue. We know that L has a zero eigenvalue with eigenvector Qp t \, 
the perturbed KMS state. This means that our task reduces to showing instability 
of all nonzero eigenvalues, and that the dimension of the nullspace of L is at most one. 

It is conventional wisdom that embedded eigenvalues are unstable under generic 
perturbations, turning into resonances. We now outline the technique we use to show 
instability of embedded eigenvalues: the PC technique. 

To do so, we concentrate first on a nonzero (isolated) eigenvalue e of Lq whose 
instability we want to show. The main idea is to construct an anti-selfajoint operator 
A, called the adjoint operator (to L), s.t. we have the following PC estimate: 

E A (L)[L,A]E A (L)>6Ei(L), (14) 

where 9 > is a srictly positive number, E A (L) denotes the spectral projector of L 
onto the interval A, and [•, •] is the commutator. Here, A is chosen to contain the 
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eigenvalue e but no other eigenvalues of Lq. Equation (14) is also called a (strict) 
Mourre estimate. If it is satisfied, then one sees that L has no eigenvalues in A by 
using the following argument by contradiction: suppose that Lip = e'tp, with e' G A 
and H^ll = 1. Then we have Ea(L)i/j = ip, and the PC estimate (14) gives on one 
hand (ip, [L, A]ip) > 9. On the other hand, formally expanding the commutator yields 

{A [L, Aty) = ty, [L - e', A]ip) = 2Re ((L - e>, Aip) = 0, (15) 

which leads to the contradiction 6 < 0, hence showing that there cannot be any 
eigenvalue of L in A. 

This formal proof is in general wrong. Indeed, both operators L and A are un- 
bounded, and one has to take great care of domain questions, including the very 
definition of the commutator [L, A]. 

Relation (15) is called the Virial Theorem, and it can be made in many concrete 
cases rigorous by approximating the hypothetical eigenfunction ip by "nice" vectors. 
The situation in which this works is quite generally given by the case where [L, A] 
is bounded relatively to L, which is in particular satisfied for iV-body Schrodinger 
systems, and systems of particles coupled to a field at zero temperature. However, in 
our case the condition is not satisfied, and as mentioned above, we have to develop a 
more general argument of this type. 

The treatment of the zero eigenvalue is similar, except that we prove (14) only on 
~RdnE^{L)P- L , where P is the rank-one projector onto the known zero eigenvector 
of L, and P 1 - is its orthogonal complement. 

2 Main results 

Our main technical result is the abstract PC estimate, Theorem 2.1. This result is 
the basis for the spectral analysis of the Liouvillian, as explained above. We point 
out that the PC estimate holds for infra-red behaviour of the form factor (see (7)) 
characterized by p > 0, which covers the physical case p — 1/2. 

Theorem 2.2 characterizes the spectrum of the Liouvillian in view of the property 
of RTE. To prove this result, we combine the PC estimate with a Virial Theorem type 
argument. It is for the latter that we need presently the more restricting infra-red 
behaviour p > 2. We think that our method can be improved. 

A direct consequence of Theorem 2.2 is Corollary 2.3 which says that the system 
exhibits RTE (recall also Proposition 1.1). 

All the results hold under assumption of the Fermi Golden Rule Condition, (18) 
and (19). In Theorem 2.4, we give explicit conditions on the operator G and the 
form factor g so that the Fermi Golden Rule Condition holds. We start by explaining 
this condition. In the language of quantum resonances, it expresses the fact that 
the bifurcation of complex eigenvalues (resonance poles) of the spectrally deformed 
Liouvillian takes place at second order in the perturbation (i.e. the lifetime of the 
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resonance is of the order A~ 2 ). 

As we have mentioned above, the Liouvillian corresponding to the particle system 
at positive temperature is given by L p = H p ® 1 — 1 <g) H p , acting on the Hilbert space 
'H p ®'Hp, so L p has discrete spectrum given by <j(L p ) = {e = Ei—Ej : E i: Ej e a(H p )}. 
For every eigenvalue e of L p , we define an operator T(e) acting on the corresponding 
eigenspace, RanP(L p = e) C H p <8> by 

T(e) = / m*(u, a)P(L p ^ e)5(L p — e + u)m(u, a), (16) 

where 5 denotes the Dirac function, and where the operator m is given by 

m(u, a) = G t gi(u, a) - G r g 2 (u, a). (17) 

Recall that g± t 2 and were defined in and before equation (10). 

It is clear from (16) that T(e) is a non-negative selfadjoint operator. The Fermi 
Golden Rule Condition is used to show instability of embedded eigenvalues. For 
nonzero eigenvalues, the condition says that T(e) is strictly positive: 

fore^O, 7 e := inf a (T(e) \ Ran P[L P = e)) > 0. (18) 

We show in Theorem 2.4 that T(0) has a simple eigenvalue at zero, the eigenvector 
being the Gibbs state of the particle system, (see (21)). This reflects the fact 
that the zero eigenvalue of L survives the perturbation, however, its degeneracy is 
removed, i.e. the zero eigenvalue of L is simple. The Fermi Golden Rule Condition 
for e = requires strict positivity on the complement of the zero eigenspace of T(0), 
i.e. 

7o := infer (r(0) \ Ran P(L P = 0)P^j > 0. (19) 

Here, P n p is the projection onto CQ P ^, and P^ v = 1 — P^p. We give in Theorem 2.4 
below explicit conditions on G and g{k) s.t. (18) and (19) hold. 

Here is our main result. 

Theorem 2.1 (Positive Commutator Estimate). Assume the IR and UV 
behaviour (7), with p > 0. Let A be an interval containing exactly one eigenvalue e 
of L and let h G be a smooth function s.t. h = 1 on A and supp/ifl o~(L p ) = {e}. 
Assume the Fermi Golden Rule Condition (18) (or (19)) holds. Let j3 > j3o, for any 
fixed < Po < oo. Then there is a \ > (depending on flo) s -t- if® < \M < ^o, then 
we have in the sense of quadratic forms on D(N 1 / 2 ) (see remark 1 below), for some 
explicitely constructed anti-selfadjoint operator A: 

h(L)[L,A]h(L) > \\^™h{L) ( 7e (l - - C(A 1/200 )) h(L). (20) 



8 



Notation. Let s be a real variable. Then 0(s) stands for a family T s of bounded 
operators depending on s, satisfying lim s „» ll^sll/ 5 — C < oo. In (20), s = A 1 / 200 . 

Remarks. 1. iV = dr(l) is the number operator in the positive temperature 
Hilbert space (see also (3) and (89)), and Pn 0o is the projector onto the span of 
fi^o, the (3-KMS state of the uncoupled system (see (4)). Also, 8 e $ is the Kronecker 
symbol, equal to one if e = and zero else. 

2. We show in Theorem A. 2 that L is essentially selfadjoint on a dense domain in 
the positive temperature Hilbert space. 

3. The commutator [L, A] is by construction in first approximation equal to N (see 
Section 7), and h(L) leaves the domain V(N 1 / 2 ) invariant (see e.g. [M]), so that (20) 
is well defined. 

4. There is no smallness condition on the interval A (apart from it only containing 
one eigenvalue of L ). 

Theorem 2.2 (Spectrum of L). Assume the IR condition p > 2 (see (7)). Let 
> Ah for any fixed 0< ( 5 <oo ;/ 3<oo. Then the Liouvillian L has the following 
spectral properties: 

1) Let e 7^ be a nonzero eigenvalue of L , and suppose that the Fermi Golden 
Rule Condition (18) holds for e. Then there is a Ao > (dependent on (3o) s.t. 
for < |A| < Ao, L has no eigenvalues in the open interval (e_,e+), where e_ 
is the biggest eigenvalue of L smaller than e, and e + is the smallest eigenvalue 
of Lq bigger than e. 

2) Assume the Fermi Golden Rule Condition (19) holds for e = 0. Then there is 
a A > (dependent on fl ) s.t. if < |A| < A and < < A , then L has 
a simple eigenvalue at zero. 

Remark. Theorem 2.2 shows that if the Fermi Golden Rule Condition holds for 
all eigenvalues of L , then L has no eigenvalues, except a simple one at zero. 

Corollary 2.3 (Return to Equilibrium). Suppose the IR condition and the 
condition on (3 as in Theorem 4-2, and that the Fermi Golden Rule Condition is sat- 
isfied for all eigenvalues of Lq. If |A| > is small (in the sense of Theorem 4-2, 
2), then every normal state w.r.t. the {3-KMS state (the zero eigenvector of L) 
exhibits return to equilibrium in an ergodic mean sense. 

The Corollary follows immediately from Theorem 2.2 and Proposition 1.1, where 
the ergodic mean convergence is defined by (13). 

Theorem 2.4 (Spectrum of T(e)). Set T p (e) := P(L P = e)T(e)P(L p = e). 
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1) Let e ^ 0. Then there is a non-negative number 5o = So(G) (independent of 
P,X) whose value is given in Appendix A. 2 (see before (97)) s.t. 

r„(e) > S inf ( \E i:j \ / dS(u,a) \g(\E i:j \,a)\ 2 J P(L P = e). 
{Eij¥*>} \ Js 2 J 

In particular, the Fermi Golden Rule Condition (18) is satisfied if the r.h.s. is 
not zero. 

2) r p (0) has an eigenvalue at zero, with the particle Gibbs state f#g as eigenvector: 

i 

where we recall that Z p ((3) was defined in (2). Moreover, if 
9o-= r inf \(tp n ,G(p r . 

{E m „<0} 

is strictly positive, then zero is a simple eigenvalue ofT p (0) with unique eigen- 
vector Vl P p and the spectrum o/T p (0) has a gap at zero: (0, 2g Z p )na(T p (0)) = 0. 
In particular, the Fermi Golden Rule Condition (19) holds. 







g— /3i? mn \ J 


r 5(E mn + u)\g\ 2 >0 

R3 



Remarks. 1. If e ^ is nondegenerate, i.e. if e = E mono for a unique pair (m ,n ), 
then (see before (97)) 5 = E„^ \(<Pn, GV„ >| 2 + £ m#mo \(<Pm,Gip mo )\ 2 . 

2. If iJ p is unbounded, then g = 0. Indeed, let m be fixed, and take n — > oo, then 

< and (</? n , G<p m ) — > 0, since y? n goes weakly to zero. Notice though that 
<7o > is only a sufficient condition for the Fermi Golden Rule Condition to hold at 
zero. 

3. For g > 0, the size of the gap, 2g Z p , is bounded away from zero uniformly in 
P > A), since 

tr 6 _ tr c~ fi^p 
lim inf = lim inf 

where E^ := E^ — E >0 (E is the smallest eigenvalue of H p ) and H p := H p — E > 
(the smallest eigenvalue of if p is zero). 



3 Review of previous results 

Proving the RTE property is one of the key problems of non-equilibrium statistical 
mechanics. Until recently, this property was proven for specially designed abstract 
models (see [BRII]). The first result for realistic systems came in the pioneering work 
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of Jaksic and Pillet [JP1,2] in 1996. 

In their work, Jaksic and Pillet prove return to equilibrium, with exponential rate 
of convergence in time, for the spin-boson system (i.e. an TV-level system coupled 
to the free massless bosonic field with N = 2; their work easily extends to general 
finite N) for sufficiently high temperatures. Their work introduces the spectral ap- 
proach to RTE. The analysis is done in the spirit of the theory of quantum resonances, 
using spectral deformation techniques, where the deformation is generated by energy- 
translation. The IR condition on the form factor is g(u) ~ uj p , uj — > 0, with p > —1/2, 
hence includes the physical case p — 1/2. However, there is a restriction on temper- 
ature: |A| < 1/(3. The spectral deformation technique imposes certain analyticity 
conditions on the form factor. 

The TV-level system coupled to the free massless bosonic field is also treated in 
[BFS4] , but the spectrum of the Liouvillian is analyzed using complex dilation instead 
of translation. RTE with exponentially fast rate in convergence in time is established 
for small coupling constant A independent of (5. Bach, Frohlich and Sigal adapt in this 
work their Renormalization Group method developed in [BFS1,2,3] to the positive 
temperature case. The IR condition is p > 0, which includes the physical case. 

In a recent work, Derezihski and Jaksic [DJ] consider the Liouvillian of the iV-level 
system interacting with the free massless bosonic field. Their analysis of the spectrum 
of the Liouvillian is based on the Feshbach method which is justified with the help of 
the Mourre Theory, applied to the reduced Liouvillian (away from the vacuum sec- 
tor). The Mourre theory in turn is based on a global positive commutator estimate 
for the reduced Liouvillian. The IR condition for instability of nonzero eigenvalues is 
p > 0, and for the lifting of the degeneracy of the zero eigenvalue, it is p > 1. 

The method for the spectral analysis of the Liouvillian we use employs the energy- 
translation generator in the Jaksic-Pillet glued positive temperature Hilbert space, as 
in [JP1,2] and [DJ]. We prove a Mourre estimate (PC estimate) for the original Liou- 
villian with a conjugate operator which is a deformation of the energy shift generator 
mentioned above. This method has been developed in the zero-temperature case in 
[BFSS] (for the dilation generator though). 

Our construction of the PC works for the IR condition p > 0, which includes 
the physical case. In order to conclude absence of eigenvalues from the PC estimate, 
the Virial Theorem is needed. So far, the systems for which the Virial Theorem was 
applied have always satisfied the condition that [L, A] is relatively bounded with re- 
spect to L, in which case a general theory has been developed, see [ABG] (for specific 
systems, see also [BFSS] for particle-field at zero temperature, [HS1] for TV-body sys- 
tems). We remark though that in [S], Skibsted extends the abstract Mourre theory 
to certain systems where [L,A] is not relatively bounded (but [[L, A], A] is). 

We develop in this work a Virial Theorem type argument in the case where the 
commutator [L, A] is not relatively L-bounded. This comes at the price that our 
estimates involve the triple commutator [[[L, A], A], A], and consequently, we need a 
restrictive IR behaviour of the form factor, namely p > 2. We think that this re- 
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striction coming from the part of the proof using the Virial Theorem (not the PC 
estimate), can be improved by a better understanding of the Virial Theorem. 

It should be pointed out that the Virial Theorem is an important tool of interest 
on its own, still currently under research, see e.g. [GG]. 

We finish this brief review by comparing our approach to that of [DJ] which, in 
the literature on the subject, is closest to ours. 

The main difference is that [D J] develop first the Mourre theory for a reduced Li- 
ouville operator, staring from a global PC estimate on the radiation sector. Using the 
Feshbach method, they show then the limiting absorption principle for the Liouvillian 
acting on the full space. It is our impression that this method is restricted to systems 
where a global PC estimate is valid, i.e. for positive temperatures, one cannot avoid 
using the generator of translations as the adjoint operator. 

The use of a different adjoint operator than the Jaksic-Pillet translation generator 
might be desirable, for instance in order to remove restrictive assumptions on the 
coupling functions. 

In our method, we modify the bare adjoint operator in such a way as to have a 
local PC estimate right from the start for the full (i.e. not for a reduced) Liouvillian. 
This method has the advantage that it works for various choices of the adjoint oper- 
ator, in fact, it was first developed (for zero temperatures) for the dilation generator 
in [BFSS]. It is true though that the use the translation generator greatly reduces the 
number of estimates to be performed, and this is the reason why we use it here. 

Let us also mention that our PC estimate is local in the spectral localization of L, 
but is of a "broad locality" in the sense that it holds in neighbourhoods of eigenvalues 
of L that need not be small, and are in particular independent of the coupling con- 
stant (the only restriction being that such neigbourhoods contain only one eigenvalue 
of the non-interacting Liouvillian). This means that we do not need two separate ar- 
guments to treat the regions "close" to (typically in a A 2 -neighbourhood) and "away" 
from the eigenvalues of L , as is often the case in Mourre theory, as well as in [DJ]. 

We do not claim that either of the two methods is better, both having, in our view, 
advantages and disadvantages. We do believe that our approach gives new insights 
and can open doors to new techniques to handle the problem of RTE and related 
spectral problems. 

4 Proof of Theorem 2.1: step 1. 

We prove in this section the PC estimate w.r.t. spectral localization in the uncoupled 
Liouvillian L , see Theorem 4.3. Step 2 consists in passing from this estimate to the 
one localized w.r.t. the full Liouvillian L and is performed in the next section. 

Our estimates are uniform in (3 > /3q (for any < (3q < oo fixed). For notational 
convenience, we set (3q = 1, see also the remark after Proposition A.l in Appendix 
A.l. 
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4.1 PC with respect to spectral localization in Lq 

We construct an operator B (see (27)) which is positive on spectral subspaces of L , 
see Theorem 4.3 (the main result of this section). 

On L 2 (R x S 2 ) and for t E R, we define the unitary transformation (jjftpj (u, a) = 

ip(u — t,a), which induces a unitary transformation U t on Fock space T = JF(L 2 (R x 
S 2 )): U t = T(U t ), i.e. for ip e J 7 , the projection onto the n-sector of U t ip is given 
by {U t ip) n («!,..., u n ) — ip n (ui — t, . . . , u n — t). Here and often in the future, we do 
not display the angular variables «i, . . . , a n in the argument of ip n . U t is a strongly 
continuous unitary one-parameter (t G R) group on T . Its anti-selfadjoint generator 
Aq, defined in the strong sense by dt\t=oUt = Aq, is Aq = —dT(d u ). The domain of 
the unbounded operator A , V(A ) = {ip e T : d t \t=oU t ip G J 7 }, is dense in J 7 , which 
simply follows from the fact that Aq is the generator of a strongly continuous group. 
From now on, we write U t = e tA °, t 6 I. The following result serves to motivate 
the definition of an operator denoted by [L, A ] (see (23) below). The proof is not 
difficult and can be found in [M]. 

Proposition 4.1 On the dense set V(L ) fl V(N), we have e~ tA °Le tA ° = L + 
tN + \I t , where I t is obtained from I by replacing the form factor g by its translate 
g l , and g l (u, a) = g(u + t,a). We obtain therefore 

d t \ t=0 e" tAo Le tAo =N + \I, (22) 

where I — G\ <S> (a*(d u gi) + a{d u g\)) — G r ® (a*(d u g 2 ) + a(d u g2)). The derivative in 
(22) is understood in the strong topology. 

On a formal level, we have d t \t=o e~ tA °Le tA ° = —A L + LA = [L,A ], which 
suggests the definition of the unbounded operator [L,A ] with domain V([L,A ]) = 
V(N) as 

[L,A }:= N + XI. (23) 

We point out that the operator [L, A ] is defined as the r.h.s. of (23), and not 
as a commutator in the sense of LA — A L. Remark that [L, A ] is positive on 
T>(N) nRan , where Q is the vacuum in T . Indeed, from Proposition A.l, it follows 
(take e.g. c = 1/4) [L, A ] > f N - 0(A 2 ), so that P£[L, A ]P£ > (3/4 - 0(A 2 )) P£. 
On the other hand, Pq[L, A ]P n = 0, so if we want to find an operator that is positive 
also on Cf2, then we need to modify A . 
For a fixed eigenvalue e e cr(L p ), define 

b{e) = e\(QR 2 IQ-QIRlQ) , (24) 
R, = ((L - e) 2 + e 2 y 1/2 . 

Here, 6 and e are positive parameters, and Q, Q are projection operators on TC defined 

as 

Q = P(L P = e)®P n , Q = l-Q. (25) 
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In what follows, we denote R e := QR e . 

Proposition 4.2 The operator b = 6(e) is bounded and [L, b] = Lb — bL is well 
defined on V and it extends to a bounded operator on the whole space. We denote 
the extended operator again by [L, b] . 

Proof. The operator b is bounded since both IQ and QI are bounded. Further- 
more, since ||L i? e || < 1 + |e|/e and ||L Q|| = |e|, then [Lq,6] is bounded. Moreover, 
since ||/Q|| < C and \\IR 2 JQ\\ < Ce- 2 \\(N +1)IQ\\ < 2Ce- 2 \\IQ\\ < Ce~ 2 , then also 
1 1 [I,b] 1 1 < oo. We used the fact that Ran/Q C Ran P(N < 1), since I is linear in 
creators and NQ = 0. ■ 

We define the operator [L, A] by V([L, A]) = V(N) and 

[L, A] := [L, A ] + [L,b]=N + \I+ [L, b] . (26) 

Again, we point out that [L, A] is to be understood as the r.h.s. of (26) (with [L, b] 
defined in Proposition 4.2). The commutator notation [L,A] is chosen because in 
the sense of quadratic forms on V(L ) n V(N) fl V(A ), one has (</?, [L, A](p) = 
2Re (Lip, Aip) with A = A + b. Define now the operator B by V(B) = V(N) and 

B := [L, A]-±N = ^N + XI + [L, b}. (27) 
Here is the main result of this section: 

Theorem 4.3 Let e G o~(L p ) and let A be an interval around e not containing 
any other eigenvalue of L p . Let be the (sharp) indicator function of A and set 

= Ea(Lq). Assume that the Fermi Golden Rule Condition (18) (or (19)) holds. 
Then there is a number s > s.t. if < 9, e, eO^ 1 , 9\ 2 e~ 3 < s, then we have on 
V(N 1 / 2 ), in the sense of quadratic forms: 

E° A BE° A > — le E A (l - l5 efi Pn fj , )E A , (28) 
where Pn [j0 is the projector onto the span ofVLp$ defined in (4). 

An essential ingredient of the proof of Theorem 4.3 is the Feshbach method, which 
we explain now. 

4.2 The Feshbach method 

The main idea of the Feshbach method is to use an isospectral correspondence between 
operators acting on a Hilbert space and operators acting on some subspace. We 
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explain this method adapted to our case. For a more general exposition, see e.g. 
[BFS2] and [DJ]. 

Consider the Hilbert spaces TL e defined by TC e = R&n Xu E A , where Xu = x(N < v) 
is a cutoff in N, and v is a positive integer. With our definitions of Q,Q, (see (25)) 
we have 

H e = Ran X uE° A Q © Ran X uE A Q. (29) 

Define Q\ = XuE A Q and Q 2 = XvE A Q an d se ^ Bij = QiBQj, i,j = 1,2. The opera- 
tors Bij are bounded due to the cutoff in N. Notice that Qi,2 are projection operators 
(i.e. Q\ 2 = <3i,2) since Xv commutes with E A and Q. 

The main ingredient of the Feshbach method is the following observation: 

Proposition 4.4 (isospectrality of the Feshbach map). If z is in the resolvent 
set of B 22 (i.e. if [B 22 — z )~ x \ RanQ2 exists as a bounded operator) and if 

\\Q 2 {B 22 - z)- 1 Q 2 BQ 1 \\ < 00, \\QiBQ 2 {B 22 - z^Qzj < 00, (30) 

then we have z e o~#(B) •<=>- z G a#(£ z ), where the Feshbach map £ z = £ Z (B) is 
defined by B 1— > £ z = Bu — B\ 2 [B 22 — z)~ 1 B 2 i, and a# stands for a or a pp (spectrum 
or pure point spectrum). 

The proof of Proposition 4.4 is given in a more general setting e.g. in [BFS2], 
[DJ], we do not repeat it here. We use the isospectrality of the Feshbach map to show 
positivity of B in the following way (see also [BFSS]): 

Corollary 4.5 Let i? = inf o~(B \ 7i e ) and suppose that B 22 > q 9Q 2 for some 
d > —00, and that inf o~(£$) > S uniformly in d for 1? < #1, where S and $ 1 are 
two fixed (finite) numbers. Then we have $0 > min-fi?, inf a(£# )}. 

Remarks. 1. All our estimates in this section will be independent of the N- 
cutoff introduced in (29). In particular, $ , S are independent of v. This will 
allow us to obtain inequality (28) on V(N l l 2 ) from the corresponding estimate on 
Ranx(iV < v) by letting v — > 00 (see (50) below). 

2. The condition inf o~{£$) > S uniformly in d for 1? < implies that $0 7^ —00. 

Proof of Corollary 4-5. If $0 > 1?, then the assertion is clearly true. If $0 < 
then i? is m the resolvent set of B 22 , and it is easy to show that (30) holds for z = $0, 
so $ e cr(^o)> i- e - ^0 > inf cr(£^ ). ■ 

4.3 Proof of Theorem 4.3 (using the Feshbach method). 

We apply Corollary 4.5 to the operator 

B' = B-5 efi 5P^ o , (31) 
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where 8 e $ is the Kronecker symbol, i.e. 5 ej o is one if e = and zero else. The positive 
number 5 will be chosen appropriately below, see after (48). 

First, we show that B 22 > (3/4 — b~ e flb~)Q2 (see (33)), then we show that £$ > 
— 1 — 5 e fiS =: S (see Proposition 4.6), uniformly in ^ for ^ < 1/2 — 5 e $5. Invoking 
Corollary 4.5 will then yield the result. Notice that due to the cutoff \v in (29), Bij, 
i,j G {1, 2} are bounded operators. All the following estimates are independent of v. 

We first calculate B' 22 = Q 2 B'Q 2 . Using QQ2 = 0, and 5 e>o Pn 0o Q 2 = o~ e ,oQ2, we 
obtain from (31) and (27) 



(32) 



B'22 = Q2 (JqN + XI + 9X 2 (R 2 JQI - IQIR]) - 5 efi 5^j Q 2 . 
Proceeding as in the proof of Proposition A.l, one shows that Vc > 0, 
(V, A/>) J < c\ \N l / 2 4>\\ 2 + C^\\d u9l \\h I M I 2 - 

With our assumptions on g, \\d u gi\\ 2 L2 < 00, uniformly in f3 > 1. Using the inequality 
above with c = 1/10 and ||P e /(5/|| < Ce~ 2 , we obtain 

B'22 > Q2 (JqN - 0(X 2 + 9X 2 e~ 2 ) - 5 efi d^j Q 2 . 

As can be easily checked, Q2 = Q2PQ, so we have NQ 2 > Q2, and we conclude that 
there is a si > s.t. if A 2 + #A 2 e~ 2 < si, then 



(33) 



(34) 



B'22 ^{jq- 6 ^ 6 ~ °( X2 + ^V 2 )) Q 2 > ^ - 5 efi 5^j Q 2 . 

In the language of Corollary 4.5, this means we can take d = 3/4 — <5 ej o^. 
In a next step, we calculate a lower bound on £# for -d < 1/2 — S ei oS. 

Proposition 4.6 We have, uniformly in d for ■§ < 1/2 — 5 ei0 S: 

£0 > 2tt^(1 - 59)Q 1 (Y(e) - ^P^ - 0(e^ + eO' 1 + £AV 3 )) Q u 

where the error term is independent of 5. Recall that Q P p is the particle Gibbs state 
defined in (21). 

Proof of Proposition 4-6 By definition, £# = B' u — B' l2 (B' 22 — d)~ l B' 21 . We show 
that B' u is positive and B' l2 (B' 22 — d)~ l B' 2l is small compared to B' n . 

With QQi = 0, QQi = Qi and deflP^^i = 5 efi P^ v Qi, we obtain from (31) and 
(27): 

B'n > 26X 2 Q 1 (lR 2 J - Qi ~ 0(X 2 ), (35) 
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where we used XI > -±N - 0(A 2 ) and QiN = 0. 

Let us now examine B[ 2 (B' 22 — i?) -1 .^. Notice that from (32), we get 

Q 2 (B' 22 - d)Q 2 = ^Q 2 N 1/2 (1 _ ip (7? + s^N' 1 + K 1 )N 1 / 2 Q 2 , (36) 

where we defined the bounded selfadjoint operator K\ acting on RanQ 2 as 

k x = y n ' 1/2 ( x! + ^(R^Q 1 - iQiR])) n~ 1/2 . (37) 

Since ||Q 2 W~ 1/2 || < 1 and | \I(N + l)- l ' 2 \ \ < C, we get ||Xi|| < C(A + tfAV 2 ). Now 
on RanP^-, we have N > 1, so since we look at s.t. + 5 et0 5 < 1/2, we obtain 

l-^ + ^A^l-y^- (38) 
Therefore we can rewrite (36) as 

Q 2 {B' 22 -$)Q 2 = ^g 2 7V 1 /2(i _m(j) + 6eo6)N -if 2 ( 1 + K2) 

x (1 - f (tf + 5 e ,o5)iV" 1 ) 1/2 iV 1/2 Q 2 , (39) 

where K 2 = (l - f(0 + c^)^- 1 )" 172 ^ (l - f (0 + S^N- 1 ) " 1/2 , and ||X 2 || < 
f H^iH = 0(A + ^A 2 e- 2 ) « 1. We have thus from (39): 

Q 2 (B' 22 - i?) -1 Q 2 = ^g 2 7V- 1 / 2 (]l-f^ + 5 e) o5)7V-i)- 1/2 X 2 

x (11 - f (0 + 5 e ,o5)iV^) _1/2 iV- 1/2 Q 2 , (40) 

where X = (1 + K 2 y xj2 is bounded and selfadjoint with ||K|| 2 = \\K 2 \\ = + 
K 2 )^ 1 \ \ < < 2. We have therefore, from (40) and (38), and uniformly in for 

< 1/2 - 5 e , 5: 

<^,^ 2 (fi 22 -^)- 1 5 21 ^> = ^||K(]l-f(^ + 5 e , 0( 5)iV-i)- 1/2 iV-V 2 ^||2 

< 2^||iV-V^|| 2 = 5||iV- 1 /2 J B 21 ^|| 2 . (41) 

Notice that I^ 2 = B 12 and _B 21 = B 21 . Now, remembering (27), and since NQ 1 = 
and Q 2 Q = = QQi, 

AT 1 ^ = jv -i/2q 2 [a/ + ^ A ( Lq _ e )R 2 j _ ^AR 2 /(L - e) + flA 2 (/i? 2 / -i^JQJ^Qx. 

Using ||JV-V2Q 2 || < 1, H/Qill < C, H/Qxll < C, WN^IW < C, (L - e)Q 1 = 0, 
\\(Lq- e)R e \\ < 1, we get 

I IN-^B^W 2 < C{\ 2 + fl 2 AV 4 ) 1 1^| | 2 + 2£ 2 A 2 \\R e IQiip\ | 2 , 
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thus with (41), we obtain 

- B' 12 (B' 22 - tfy'B'^) > -109 2 \ 2 (ip, Q 1 IR 2 e IQ 1 ^ - 0(A 2 + # 2 A 4 e" 4 
and so, together with (35), we get, uniformly in for i9 < 1/2 — 8 e:0 S: 

E* > 26\ 2 (l - 59)Q 1 (IR 2 e I - ^PfyQi ~ 0(\ 2 + £ 2 AV 4 ). (42) 

We point out that the error term in the last inequality does not depend on S. With the 
choice of parameters we will make (see (68)), (42) shows that £# > — 1— 5 ei0 5 uniformly 
in $ for d < 1/2 — S ej0 S, i.e. in the language of Corollary 4.5, S = — 1 — 5 e ,oS- 

The remaining part of the proof consists in relating the strict positivity of the 

2 

nonnegative operator Q 1 IR e IQi to the Fermi Golden Rule Condition. We let I a 
and I c = I* denote the parts of / containing annihilators and creators only, so that 
I = I a + I c . Thus 

Q 1 IR 2 JQ 1 = QjXhQi = QJ a R%Qi. (43) 

In the first step, we used I a Qi = and Q±I C = (since I a Pn = 0) and in the second 
step, we used Q\I a Q = Qih (since I a Q = 0). Now write 



QihRJcQi = Qi m*(u,a)a(u,a)R 2 e (e)a*(u,a')m(u',a') Q u (44) 



where m is defined (17), and where we display the dependence of R\ on e. The 
operator- valued distributions (a and a*) satisfy the canonical commutation relations 
[a(u, a), a*(u f , a')} = 5(u — u')5(a — a'). Next, we notice that the pull-through formula 
a(u, a)Lf — (Lf + u)a(u, a) implies 

a(u, a)R 2 (e) = R 2 (e — u)a{u, a). (45) 

Using the CCR and formula (45) together with the fact that a(u, a)Q± = 0, we 
commute a(u, a) in (44) to the right and arrive at 

(44) = Qi J m*(u, a)R 2 (e - u)m(u, a) Q x . (46) 

We can pull a factor Pq out of Q\ and place it inside the integral next to R 2 (e — u) and 
thus replace R 2 (e — u) by ((L p — e + u) 2 + e 2 )^ 1 . Notice that e((L p — e + u) 2 + e 2 y l — > 
5(L P — e + u) as e — > 0. More precisely, we have 

Proposition 4.7 There is an s 2 > s.t. for < e < s 2 , we have 

Qi J m*(u, a) ((L p -e + u) 2 + e 2 )' 1 m(u, a) Q 1 > Q^- (r(e) - 0{e x ^)) Q x . 

18 




Proposition 4.7, which we prove in Appendix A. 3, together with (42)-(44) and (46) 
yields (34), proving Proposition 4.6. ■ 

Now we finish the proof of Theorem 4.3. If the Fermi Golden Rule Condition (18) 
holds, then for e ^ 0, we have T(e) > 7 e > on RanQi, so we obtain from (34), and 
under the conditions on the parameters stated in Theorem 4.3: > 7r^7 e , so by 
Corollary 4.5: 

inf a(B \ H e ) > min{l/2, nOXh' 1 ^} = vr— 7e , (47) 

since by our choice of the parameters (see (68)), we will have ^ < (27r 7e ) _1 . 

For e = 0, we have T(0) = r(0)P£ P , since T(0)^ = (see Theorem 2.4), so 

Proposition 4.6 gives 



({70 - ^} Pfy - 0(e^ + dT 1 + 9X 2 e^ Ql . (48) 



For some fixed < a < 2 (2-i) (independent of 9, A,e), there is a S3 > s.t. if 
< 9\ 2 e~ 1 < s 3 , then 70 — ^ff? > —a, which gives with (48): 

9 \ 2 / \ 
£4 > 7T Ql (-aP- $ -0(e^ + ee-i + e\>e-*)) Ql 

> 7r ^(_ a _o( e 1 /4 + er i + ^ A V 3 ))g 1 

9\ 2 „ 

> —2tt aQi. 

e 

The last step is true provided e 1 / 4 + e9~ x + 9\ 2 e~ 3 < S4, for some small S4 > 0. 
Remembering that B' = B — 5Pq /}0 , we obtain from Corollary 4.5 



inf cr ((B - 5P£ 0Q ) \ Ho) > min{l/2, -2na9\ 2 /e} = -2na 



9X' 



from which we conclude that if the condition on the parameters given in Theorem 4.3 
is satisfied with s = min(si, S2, S3, S4), then 



(9\ 2 \ 
-27ra— + 5P^)E° A Xu 



= 2^ 7 oxX (1 - a(n - l)/ 7o - (1 + a/ lG )P Qp o ) E° A Xv 
9X 2 

> — 1oXu E a (1-IP^ )E° aX u, (49) 
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where we used a/70 < 2 ( J-i) ■ Estimates (47) and (49) yield V0: 

(V>, XuE° A BE° A Xv^) > — 7e XX(1 " ¥e, O Pn , o )E o A x^) ■ (50) 

Suppose now if> e V^N 1 ' 2 ). Then, since (N + l) _1 / 2 .B(iV + l)" 1 / 2 is bounded (see 
the definition of B, (27)), and since Xu — ¥ 1 strongly as z/ — > 00, we conclude that 
VV> G ©(iV 1 / 2 ): 

<^, E° A BE° A i>) > — 7e (V>, - |o e , P^, )^> , 
which proves Theorem 4.3. ■ 



5 Proof of Theorem 2.1: step 2. 

We pass from the positive commutator estimate w.r.t. L given in Theorem 4.3 to one 
w.r.t. the full Liouvillian L, hence proving Theorem 2.1. The essential ingredient of 
this procedure is the IMS localization formula, which we apply to a partition of unity 
w.r.t. N. Then, we carry out the estimates on each piece of the partition separately. 



5.1 PC with respect to spectral localization in L 

Let 1 = xii x ) + xl( x )i x e Xi e QTQO' !])> be a C°°-partition of unity. For some 
scaling parameter <j » 1, define x« — Xi(N) = Xi{N/<j), % — 1,2. The reason why 
we introduce the partition of unity is that Ixi = O (a 1 / 2 ) is bounded. Since the x% 
leave ©(iV 1 / 2 ) invariant, then [xi, [xi,B]\ = x 2 B — 2xiBxi + Bx 2 is well defined on 
V(N 1 / 2 ) in the sense of quadratic forms, and by summing over i — 1,2, we get the 
so-called IMS localization formula (see also [CFKS]): 

B = Y,XiBXi + \\Ki,\x.i,B\\. (51) 

1,2 

Furthermore, we obtain from (51) and (27), in the sense of quadratic forms on 
V(N 1 / 2 ): 

h(L)[L,A]h(L) = ^h(L)Nh(L) + J2 h ( L ^ B ^ML) + lh(L)[ Xl ,[^B]]h(L). (52) 

iU 1,2 1 

In Propositions 5.1-5.3 below, we estimate the different terms on the r.h.s. of (52). 
Then we complete the proof of Theorem 2.1 by choosing suitable relations among the 
parameters 6, A, e, a (see (68)). 
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Proposition 5.1. There is a s 5 > s.t. if X 2 a 1 < s 5 , then 



hXiBxih > ^hx 



lh. 



(53) 



Proo/. Recall that B = + XI + [L, b}. S ince Qx2 = and QIX2 = (see also 
end of proof of Proposition 4.2), we have V0: {^, X2 [L,b} X2 ^) = 0. Furthermore, 
Proposition 6.1 gives Vc > 0, XI > cN — 0(A 2 /c), so 

i>, x 2 (9iV/10 + XI)X2^) > (A X2 [(9/10 - c)N - 0(A 2 /c)] X 2^> > \a &l>) , 

provided X 2 a < s 5 and where we picked the value c = 1/10 and used X2NX2 > ^xl-' 

Proposition 5.2. We /iawe 



1 /Q \ 2 r a\2 

h Xl B Xl h+-hNh > — le (l-O(Xa 1 / 2 ))hx 2 1 h--— l0 6 efi hP^ h 

Q\2 

{e6- l + ea l / 2 + Xae- l )h 2 . 



Proof. Let F£, := F A ,(L ), where A' is an interval whose interior contains the 
closure of A, and F A > is a smooth characteristic function with support in A', s.t. 
E a (Lq)F^, = 0, where we denoted 1 — F A , =: F^,. We take A' to contain only one 
eigenvalue of cr(L ), namely e, so that (28) in Theorem 4.3 holds, with E A replaced 
by E° A ,. We have 



h X iB X ih+—hNh 



h Xl F A ,BF A , Xl h (54) 
+ ^-hNh + h Xl F^BFf /X ih+ adjoint (55) 



+h Xl F^BF° A , X ih. 



(56) 



First, we show that (55) and (56) are bounded below by small terms. To treat (55), 
notice that 



Xl F° A ,BF A , X i = X iF° A ,(9N/lO + XI+[L,b])F A , X i 



= ^ X lF A/ FlN + Xl F AI (XI+[L,b])F° A , Xl 

> Xl F° A ,(XI+[L,b])Fj, X i. 
Now for 12 G U^N 1 / 2 ), we have for any c > (see Proposition A.l) 



(57) 



01, A/02 



< A 



01,402 



+ 



(f>2,Ia(f>l 

< cxdi^w ||7v 1/2 2 || + ||0 2 ||||^ 1/2 0i||) 

< CAV 1 (||0 1 || 2 + ||0 2 || 2 )+ C (||iV 1 /V 1 || 2 + ||iV 1/2 02|| 2 ). 
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With 0! = F A , X iip, 02 = F A , X iip, this yields Vc > 0: 

>, Xl F°, A/F^) | < C7^2| | 2 + 2c\ |7V 1/2 X iVl | 2 , 

so XiF£,A/fJxi + adjoint > -4 (c^x? + cNj . Taking c < 1/40 gives then 
1 
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hNh + hxiF° A ,\IF%xih + adjoint > (1/10 - 4c)/iiV/i - CA 2 /ix^ 



> -CX 2 h X \h. 



(58) 



Next, using QFjy = and (L — e)Q = 0, we calculate 



XiF%[L,b]F% X i = X iF A ,[L -e,b}F AlXl + X Xl F^[I,b}F AlXl 
= 9\xiF A ,QIR 2 e (L -e)Ff,xi 



+9X 2 Xl F° A , \-RJQI - IQIR e + QIRJ) F% X i 



(59) 



where we used ||i? e F^,|| < |A'| 1 < C and \\I X i\\ < Ca 1 / 2 . Next, since supp/i fl 



suppF^, = 0, then Xl F A ,h(L) = Xl F A ,(h(L) — h(L )), so by using the operator 
calculus introduced in Appendix A. 4, we obtain 



Xl F A ,h(L) = X iJ dF A ,(z)(L - zy'XIiL - z )" 1 / l (L) = OtXa 1 ' 2 ) 



(60) 



From (59), we then have hxiF A ,[L,b]F A ,xih > -C^ea 1 ' 2 + Xae^h 2 , which, to- 
gether with (58) and (57) yields 



(55) > -C U -^{e9- 1 + ea 1 ' 2 + Xae^h 2 . 



(61) 



Our next step is estimating (56). Again, using QF A , = 0, we get 

Xl Ff,BFf, Xl = Xl ^(9^/10 + A/)^xi-^A 2 xi^(^/g/ + /gffi £ 2 )^;Xi 
> -C(A 2 + #A 2 ), 



where we used XI > -cN - 0(X 2 /c) and \\F A ,R 2 \\ < |A'|~ 2 < C. We thus obtain, 
since « 1: 

n\2 

(56) = h X iF%BF%xih > ~C—-h 2 . (62) 

e u 
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Finally, we investigate the positive term (54). By sandwiching (28) in Theorem 4.3 
(with E° A replaced by E\,) with Fjy, and noticing that F^,E\, = Fjy, we arrive at 

Q\2 

hxiFlBF^x.h > <K— le h Xl Ft, (l - f ^o^Vo) F% X ih 

a\2 



= 6 -^leh (xl (l - - |<W) h 



> —leh (xl (l - 2ift) - §^Pn M ) h 

g\2 

> —lehixKl-CXa^-p^PnJh, (63) 

where we used (60) in the last step once again, and — 2x\ (F^i) 2 Pn > —2Pq /30 in 
the second step. Combining (63) with (61) and (62) yields Proposition 5.2. ■ 

Proposition 5.3. We have J2i,2 h iXi, [Xi,B]]h = ^ O^"^ 1 ^ 3 / 2 )/* 2 . 

Proof. Notice that X i an d 1 ~~ X2 have compact supports contained in [0,2]. 
Now in the double commutator, we can replace X2 by 1 — X2 without changing its 
value. So if suffices to estimate [ X (N/a),[ X (N/a), B]], where X e Co°([0,2]). We 
have [x(N/a), [ X (N/a),B}} = [x(N/a), [ X (N/a), XI + [L,b}}}. It is not difficult to see 
that we have in the sense of operators on V(N 1 / 2 ): 

[X{N/<J)MN/<J)XI\\ = dx(z) J dxiOWa-zr^N/a-C)- 1 

xI(N/a - z)-\N/a -CY 1 - (64) 

We used the operator calculus introduced in Appendix A. 4. Now since \\I(N/o~ — 
z)~ 1/2 \\ < C\\(N + lf/ 2 (N/a - ^) — 1 1| < Ca^llmz]- 1 , which follows from 

sup , V ^ TI < Ca^llmzl' 1 , 
x >o \x/a — z\ 

we conclude that 

[ X (N/a),[x(N/a),\I]} <C\a^l 2 . (65) 
Next, write for simplicity X instead of x(iV/cr), and look at 

[ X , [x, [L, &]]] = 9X[x, [x, [L,R 2 JQ]}} + adjoint. 

We claim that 

[ X , [LMllQW = 0. (66) 
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Write first [L,R 2 JQ] =R 2 t [L ,I]Q+X[I,R 2 JQ]. Then [ X ,R 2 [L , 1]Q] = [%r][L q , I]Q] = 

2 2 

xR t [Lo, I]Q — R t [L , 1]Qx- Here, x — 1 ~~ X- Notice that Qx = 0, and since 
RanR 2 e [L , 1]Q C RanP(iV = 1), we have also %R^[L ,/]Q = 0, for a > 2. Simi- 
larly, [ X , [I,R 2 JQ}} = 0, so (66) follows. 

We obtain thus from (65): [x, [x, B]\ = 0(Xa~ 3 ^ 2 ), which proves the proposition. 

■ 

Now we finish the proof of Theorem 2.1. The IMS localization formula (52) to- 
gether with Propositions 5.1, 5.2, 5.3 yields 

n \ 2 cr n \ 2 

h[L,A]h > — le {l-O{\a l ' 2 ))hx 2 1 h + -hx 2 2 h--— 1 ,5 efi hP n[)0 h 
9X 2 



e -O (eO- 1 + ea 1 ' 2 + Xae' 1 + e9~ x X^a^' 2 ) h 2 . 

The sum of the first two terms on the r.h.s. is bounded below by (l — 0{Xa 1 / 2 )) h 2 , 
so we get 

OX 2 r 

h[L,A]h > —h^il-IS^Pn^-OiXa 1 / 2 )) 

-O (eO" 1 + ta 1 ' 2 + Xae- 1 + e9~ X A~ V" 3 / 2 ) ] h. (67) 

Finally, we choose our parameters. Let e = A e / 100 , a = A~ CT / 100 , 9 = A 61 / 100 , and choose 

(e,<7,0) = (44,55,26). (68) 

It is then easily verified that for small A, the conditions on the parameters given in 
Theorem 4.3 and Proposition 5.1 hold, and furthermore, (67) becomes 

h[L, A]h > X 182/100 h [ 7e (1 - §* e ,oPn M - 0(A 145 / 200 )) - 0(A 1/200 )] h 

> X^h(^(l-55 efi Pn Pfi )-O(X^ 200 ))h. M 



6 Proof of Theorem 2.2 

We follow the idea of the Virial Theorem, as explained in Subsection 1.3: Assume ip 
is a normalized eigenvector of L with eigenvalue e. If e = 0, we assume in addition 
that ip G RanPj^. Let a > and set f a := a' 1 f(iaA ), where / is a bounded 
C°°-function, such that the derivative /' is positive and s.t. f'(0) = 1 (take e.g. 
/ = Arctan). Set 

f' a := f(iaA ), and h a := y/ft. 

Furthermore, set f" := f"(iaA ). For v > and g G C£°(— 1, 1), define tp v = gi^uN)^. 
Here, a, v will be chosen small in an appropriate way. We define the regularized 
eigenfunction ip a ^ = h a ip u . Notice that 

i>a,u — > as a, i/ -> 0. (69) 
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Set for notational convenience in this section 

K := [L, Aq] = N + XI. 

The strategy is to show that (K)^ a ^ := (tp atU , Kip a>v ) — > 0, as a, v — > (see next 
Subsection, (74)). For this estimate, we need the restrictive IR behaviour p > 2, see 
after Proposition 6.1. Using the PC estimate, Theorem 2.1, we also show that (K)^ v 
is strictly positive (as a, v — > 0), see Subsection 6.2, (86). The combination of these 
two estimates yields a contradiction, hence showing that the eigenfunction ip of L we 
started off with cannot exist. 

In the case e = 0, we need to use that the product Pn go Pa l3X is small, which is 
satisfied provided j3\\\ < C, see (11). 



6.1 Upper bound on (K)^^ 

Using (L - e)ip = and that [N, I] is iV^-bounded, we find that 

(f a (L - e)>^ = (gJ a (L - e)g u )^ = (Jag v [XI,g v % = 0(Aq,-V/ 2 ). (70) 
Next, observe that 

2Im(/ Q (L-e)>^ = ([L, if a ])^ = Re <[L, if a ])^ = Re (f a N + X[I, if a ])^ , (71) 
where we used in the last step 

[L ,i/ Q ] = J df{z)(iaA Q - z)~ 1 [LQ,A Q \(iaA G - z)' 1 
= J df{z){iaA - z)~ 2 N 

= m 

since A and N commute (second step) and we made use of (113) with p — 1 in the 
last step. The commutator [I, if a ] is examined in 

Proposition 6.1. The following equality holds in the sense of operators on V{N 1 / 2 ) 
or in the sense of quadratic forms on V^N 1 ^): 

ifa] = f a ad\ o (/) - \aZad\ (/) + R, (72) 

where we assume that the k-fold commutator ad\ o (I) := [•••[/, Aq], Aq, •• • , Aq] is 
N 1 / 2 -bounded (or N 1 ^ i -form bounded) for k = 1,2,3. The term R satisfies the esti- 
mate RN- 1 ' 2 , N-^RN- 1 ^ = 0(a 2 ). 
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Proof. Using the operator calculus introduced in Appendix A. 4, we write 

[I,if a ] = J df(z)(iaA - z)- l [I,A ](iaA - z)- 1 

= fLad-Aoi 1 ) ~ ia J df(z)(iaA - z)~- 2 ad 2 Ao (I)(iaA - z)' 1 

= f a ad\M) - l®f'«ad 2 Ao (I) -a 2 J df(z)(iaA - z)- 3 ad 3 Ao (I)(iaA - z)~\ 

The last integral is defined to be R, and the estimates follow by noticing that A and 
N commute. ■ 
Notice that it is here that we need ||ad^ (/)A" 1//2 || < C, k — 2,3, hence the more 
restrictive IR behaviour p > 2. We obtain from (72) and recalling that / = [I, A ]: 

(71) = Re (f' a K)^ - ^Re J))^ + 0{Xa 2 v^ 2 ) 

= (K)^ + XRe(h a [h a ,XI]- l -af:ad 2 Ao (I)^ + O^Xa 2 ^ 1 ' 2 ) 
= (K)^+0(Xa 2 u- 1 / 2 ). (73) 
We used in the last step that the real part in the second term above is 

[ha, [h a J]} - l -aX,ad 2 Ao (I)}\ = Oia'u- 1 / 2 ), 
since ad A (I) is A rl / 2 -bounded. Combining (73) and (70), we obtain 

/ Z/l/2 2 \ 

^<C\( — + - U - 2 Ml 2 - (74) 



6.2 Lower bound on (K) 



Let A be an interval containing exactly one eigenvalue, e, of L p . We introduce two 
partitions of unity. The first one is given by 

Xa + Xa = 1, 

where xa € C°°(A), x&{ e ) = 1- We localize in L, i.e. we set xa = Xa{L). The 
second partition of unity is given by 

x 2 + t = h 

where x ^ C°° is a "smooth Heaviside function", i.e. x( x ) = if a: < and x( x ) = 1 
if x > 1. We set for n > 0: Xn = x(N/n), x 2 n = \ — xi- We will choose n < 1/V, so 
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that XnVv = Xntp- The last equation will be used freely in what follows. We are going 
to use the IMS localization formula (51) with respect to both partitions of unity, and 
we start with the one localizing in N: 

( K )i, a , v = \XnKXn + Xn K Xn + ^[Xn,[Xn,K]] + ^[x n ,[X n , K ]]) 

\ ' 4>a,v 

> (K) Xn ^ + ^\\x n ^A 2 -0(Xn-^), (75) 

where we used that K > n/2 on RanP^ - , and the estimate (65) with a replaced by n. 
Next, from the IMS localization formula for the partition of unity w.r.t. L, we have 

( K ) Xn f a „ = (xaKxa + XaKxa + R) X71 ^ 

> ( X a(K + [L, b}) X A + XaKxa + R) Xn ^ v - X 19/50 O(an + Xn^ 2 ) 

> 0\\XAXn*Pa,u\\ 2 ~ CeS^oWPn^XAXn^aAl 2 + (Xa K Xa + R ) X n^ 

-X 19/50 O(an + Xn' 1/2 ). (76) 

Here, several remarks are in order. First, we have set 2R = [x A , [xa> -^]] + [xa, [xa, K]], 
and we have used in the second step the fact that 

& ])xax^, = (1 L ~ e ' b ]) X Ax n h a ^ = 2Re (Xa(L - e)h aX nip, bxAXnKf) 
= X^Oian + Xn- 1 / 2 ). 

We recall that b is a bounded operator (see Proposition 4.2), with ||6|| = 0(A 19 / 50 ). 
In the last step in (76), we used the positive commutator estimate, Theorem 2.1, in 
the following way. For e^O, Theorem 4.1. gives right away Xa{K + [L, b])xA > #Xa> 
where we recall that [L,A] = [L,A ] + [L,b], and b is defined in (24). We have set 
6 = CX 2 . In the zero eigenvalue case, e = 0, we have 

^91/50 

^91/50 5A 91 / 50 

> -^—lo\\XAXn^a,v\\ 2 g ^\\Pn^ XAXn^a,u\\ 2 ■ 

Setting again 9 = CA 91/50 yields (76). 

We now estimate the remainder term R. Notice that the same observation as 
at the beginning of the proof of Proposition 5.3 shows that we have the estimate 
( R ) X ni> a ,v = 2ilm (XAXn^ a ,v, [Xa, K \Xntp a ,v) ■ Therefore, 

( R ) Xn ^\ < CWXaXuKiPW \\[xA,K]xnh a ^\. (77) 

Now we have on V(N): [x A , K] = J dx A (z)(L — z) -1 ^, L](L — z)' 1 , where we recall 
that (L — z)~ x leaves V{N) invariant. Furthermore, 

[K, L] = X[N, I] + A[J, L ] + A 2 [J, /] = X[N, I] + XI{ud u g) + A 2 [J, J], (78) 



<^+M]> XAXn ^ > -T-<7o(l-5P^ )-0(A^oo ))? 
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where I(ud u g) is obtained from / by replacing the form factor g by ud u g. The last 
commutator in (78) is bounded, and the other two are iV^-bounded, so we obtain 

\\[x A ,K] X nh a ^\\ = 0{\n l l 2 )\\xni>«A- (79) 
Next, we estimate the first term on the r.h.s. of (77): 

HXaXA^II = \\{L- zY 1 Xa{L - e)Xnh a i/i\\ 

< C\\{L-e) X nKn 

< CWn-'XiN, I}XnKn + 0(\n- 3 / 2 ) + C\\ Xn (L - e)M>|| 

< C\n- 1/2 \\x'Jj a A + 0(Xn~ 3/2 + an). (80) 

Combining this with (79) and (77), we arrive at the estimate 

{R) Xn ^\ < C\ 2 \\ X ' n ^A WXn^A +0(\ 2 n- 1 + Xan 3 ' 2 ). (81) 

There is one more term in (76) we have to estimate: (XA^XA) Xn f a „■ Since 
P^{N + XI)P£ > and since P n IPn = 0, we have the bound K > P^XIPq + adj. > 
— CX, which implies 

(XaKXa)^ > -CX\\x A Xn^ a A 2 - (82) 
Using (82) and (81), we obtain from (76) 

( K ) X nf a ,v ^ 0\\Xn^a,u\\ 2 + CX)\\x A Xn^a,u\\ 2 ~ CdS e,0 1| ^, XAXn^ II 2 

-CX 2 \\ X ' n iPaJ WXn^aJ ~ X 19 / 50 O(an + An" 1 / 2 ) 

-XO(an 3/2 + Xn- 1 ). (83) 

Next, we have for any i],e > 0: 

||XnVV/|| llXnVvll ^ V\\Xn^a,u\\ 2 + V~ 1 \\x' n ' l l J a,u\\ 2 

< (67T 1 + V)\\Xn^aA 2 + V^^llXn^aA 2 - 

In the second step, we used the standard fact that we can choose the partition of 
unity s.t. HXnV'll 2 ^ e IIXnV'l| 2 + ^llXn^ll 2 ) f° r an y e > 0. Combining this with (83), 
we obtain from (75): 

> (e-cx'ier' + vmxn^aA' + in/^-cx^h-^iix^A 2 

-C9S e:O \\Pn , o XAXn^aA 2 - (0 + C^IIXAXn^H 2 

-0{Xan 3 ' 2 + A 19 / 50 an + A 69/50 n~ 1/2 ). 



Consider A small and fixed. Then if 

n 
2 



- - C^e- 2 > 9, (84) 
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we obtain 



(K) hai ,„ > 0\\h a ip v \\ 2 - C95 efi \\P npo XAXnh a iJ u \\ 2 ~ Oier]- 1 + 77 + an 3/2 + n~ 1/2 ) 
-C6(n- 1 + n~ 3 + aV). (85) 

On the last line, we used (80). Let us choose the parameters as follows: 

e = a 1 / 10 , V = a 1 / 20 , n = oT 1 ' 2 , 

then (84) is verified, and furthermore, (85) reduces to 

> V\\^A 2 ~ C05 e , o \\Pn^XAXnipaA 2 - 0(a 1/2 °). (86) 

On the other hand, recalling (74), we obtain by choosing the parameters v and a as 
v = a 3 : 

( K )^„ < C * 1/2 - (87) 
Since HV^II ~ * ll^ll = 1 as a, z/ ^ 0, and since 

-C^ e , ||P^, XAXn^,,|| 2 - -C^ e ,o||P^, P4 A VH 2 

(recall that ^ = Pq x ip if e = 0), we obtain thus for small a from (86) and (87) the 
inequality 

\{l-C8 e ^P^ gx n 2 )<Ca^. (88) 

For e ^ 0, this is a contradiction, and it shows that there can not be any eigenvalues 
of L in the interval A. Remark that there is no smallness condition on the size of 
A, except that it must not contain more than one eigenvalue of L , so we can choose 
A = (e_,e+). 

Let us look now at the case e = 0. Again, we reach a contradiction from (88), 
provided \\Pn -^% x ^\\ 2 ^ ^ ^ n ^ ms case ! we conclude that zero is a simple eigen- 
value of L. Now the fact that HPn^, P^g J| — 0(/3\X\) follows immediately from (11), 
so taking f3\X\ small enough finishes the proof of Theorem 2.2. ■ 



Acknowledgements 

The author thanks I.M. Sigal for his support and advice. Many thanks also go 
to J. Frohlich and R. Froese for stimulating discussions and to the referee for helpful 
comments. During the writing up of this work, the author has been supported by 
an NSERC PDF (Natural Sciences and Engineering Council of Canada Postdoctoral 
Fellowship), which is gratefully acknowledged. 



29 



A Appendix 

A.l Selfadjointness of L and some relative bounds 

We introduce the positive operator A = dr(|u|) with domain T>(A) = {ip G T~i : 
||A^|| < 00} and the number operator 

N = dr(l) (89) 

with natural domain V(N) = {ip G H : \\Nip\\ < 00}. 

Proposition A.l (Relative Bounds). Set L? = L 2 (RxS 2 ), and let < f3 < 00 

be a fixed number. 

1) Iff EL 2 , then\\a{f)N- l l 2 \\<\\f\\ L ,. 

2) If\u\- l ' 2 f G L 2 , then \\a(f)A-^ 2 \\ < \\ \u\-^ 2 f\\ L 2. 

3) For ip G V(N l l 2 ) and ip G ©(A 1 / 2 ) respectively, we have the following bounds, 
uniformly in f3 > (3q: 

wm 2 < ciiGiKiiiv^f + ii^ii 2 ), 

||/^|| 2 < C7||G||(||A 1 /V|| 2 + |H| 2 ). 
Here, C < C'(l + Pq 1 ), where C is independent of (3,(3 . 

4) For ip G V(N l l 2 ), any c > 0, and uniformly in (3 > (3 , we have 

1 fiA 2 f 
|(^A^)|<c||iV 1 / 2 ^|| 2 + ||G|| 2 |M| 2 / (l + ^u-'M^k. 

5) For ip G ©(A 1 / 2 ), any c > 0, and uniformly in (3 > (3q, we have 

t9\ 2 r \n\ 2 

|(^A7^)|< C ||AVV|| 2 + ^||G|| 2 ||^|| 2 / (l + /3-W)^^. 

Remarks. 1. The parameter (3® gives the highest temperature, To = l//?o, s-t. our 
estimates 3)-5) are valid uniformly in T < To. To can be fixed at any arbitrary large 
value. Since we are not interested in the large temperature limit T — > 00, we set from 
now on for notational convenience T = 1. 

2. Notice that 4) and 5) tell us that Vc > (with the O-notation introduced after 
Theorem 2.1), 

|A/| <cN + 0(A 2 /c), |A/| < cA + 0(A 2 /c), 

where we understand these inequalities holding in a sense of quadratic forms on 
V(N 1 ' 2 ) and T>(A 1/2 ) respectively. 
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Proof of Proposition A.l. The proof is standard (see e.g. [BFS4], [JP1,2]), we 
only present the proof of 3), as an example of how to keep track of f3. 

From ||/^|| 2 < 4||G|| 2 (||a*(^|| 2 +||a*(^f+||a(^l| 2 + ll«(^|| 2 ), and 
using the CCR [a*(f),a(g)\ = (f,g), we get 

IK(#i, 2 )^|| 2 = (^,0(^1,2)0*^1,2)^) = \\a(gi,2)ip\\ 2 + Ibi^ll^UV'll 2 , 

so ||/^|| 2 < 8||G'|| 2 (||a(( 7l ^|| 2 +||a(( 72 )|| 2 + 2||( ?1 || 2 2 ||^|| 2 ), where we used \\ 9l \\ L , = 
||(?2||l 2 , since gi(u, a) = —g2(—u,a). Using 1) and 2) above, we get 

||/^|| 2 < 16||G|| 2 |b 1 ||| 2 (||7V 1 />|| 2 + ||VH 2 ), 

ii/^ii 2 < leiiGfii^r^n^diA^f+ii^ii 2 ). 

Next, we show that ||<7i||i,2 < C and || |m| _1 / 2 (? 1 || L 2 < C, uniformly in (3 > (3 . Indeed, 
notice that HoiH 2 ^ = /^(l + 2fi)\g(u,a)\ 2 dudS(a) = \\g2\\ 2 L 2, where we represented g 
in the integral in spherical coordinates. Since we have 1 + 2//= 1 + 2(e l3uJ — 1) _1 < 
1 + 2f3~ 1 uj- 1 < 1 + 2Pq 1 u~ 1 , uniformly in (3 > fi , we get with (7) (for p > 0) the 
following uniform bound in (3 > (3 : 

\\9i,2\\h< 2 [ (l + Po 1 ^ 1 )\9(k)\ 2 d 3 k = C<oc. (90) 

Similarly, || — 1 ^ 2 ^i < 2 / R3 (l + f3Q 1 ou~ 1 )u~ 1 \g(u}, a)\ 2 d 3 k = C < 00, uniformly in 
(3 > (3q. It is clear from the last two estimates that C satisfies the bound indicated 
in the proposition. ■ 

These relative bounds and Nelson's commutator theorem yield essential selfad- 
jointness of the Liouvillian: 

Theorem A. 2 (Selfadjointness of the Liouvillian). Since H p is bounded be- 
low, there is a C > s.t. H p > —C. Suppose that [G, H P ](H P + C) -1 / 2 is bounded in 
the sense that the quadratic form ip 1— > 2ilm (Gip, H p ip) , defined on V(H P ), is repre- 
sented by an operator denoted [G,H p ] Q , s.t. [G, H p ] Q (H p + C) -1 / 2 is bounded. Then 
VA G M., L is essentially selfadjoint on 

V := V(H p ) <g> V(H P ) <g> P(A) CH P ®H P ® F(L 2 (R x S 2 )). 

Proof. The proof uses Nelson's commutator theorem (see [RS], Theorem X.37). 
Let N = (H p + C) <g> l p + l p <g> (H p + C) + A + 1, then N is selfadjoint on V and 
M > 1. Also, L is defined and symmetric on V . 

According to Nelson's commutator theorem, in order to prove Theorem 1.2, we 
have to show that Vip G V and some constant d > 0, 

HLVII < (91) 
\(Li{>,tfif>)-(tfi/>,Li{>)\ < d\\N 1/2 ij\\ 2 . (92) 
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Estimate (91) easily follows from HV^II < 1, WLf^W < 1 and {{IM^W < 
\\I(A + iy 1/2 \\ ||(A+ 1) 1/2 (A + l)" 1 !! < d (by 3) of Proposition 6.1). 

To show (92), notice that L commutes with M, so the l.h.s. of (92) reduces to 

| (7V, Mf) - (Afy, | < | (/^, A^) - (A^, m | + K, (93) 

where 

X = \{Ii/>,((H p + C)®l+l®(H p + C))i/>) 

- (((H p + C) <g> 11 + 11 <g> (if, + C))^, J^)| • (94) 

Let us examine the first term on the r.h.s. of (93). It is easily shown that since 
\u\9i,2 e L 2 (R x S 2 ), then a*(#i )2 )A = Aa*(g li2 ) + a*(\u\g lj2 ) on V(A). This shows 
that a#(<?i ;2 ) leave P(A) invariant and so we have Wip G £>o : 

|</^A^)-(A^,/V)I 
= |(^(/A-A/)^)| 

= |(V>, (d ® (a*(|u|^i) - a(|«|^i)) - G r <g> (a*(|«|^ 2 ) - a(|u|^ 2 )))^>| 

< cii^iiKA+i^vii^cii^Vll 2 , 

where we used Proposition 6.1 in the third step. 

Now we look at K given in (94). Using the specific form of I (see (9)), we can 
write K < \K±\ + \K 2 \, where 

K x = (G l ®(a(g 1 ) + a*(g 1 ))ij,(H p + C)®lxlj) 

- ((H p + C)® IV, Gi ® (a( 9l ) + a*( 9l ))il)) , 
K 2 = (G r ®(a(g 2 ) + a*(g 2 ))tl>,l®(H p + C)tl>) 

- (1 ® (H p + C)i/>, G r ® (a(g 2 ) + a*(g 2 ))i/>) • 

We examine i^i. Let ^ G X>o, then (if p + C) 1 ^ 2 ^ G 7i, and so 

K x = 2iIm{G l ®(a(g 1 ) + a*{g 1 )),(H p + C)®lil)) 
= 2Tm ((a( 9l ) + a*( 9l ))^, [G, H P U) , 

so we obtain \K X \ < c||(A + \\(H p + Cf' 2 ® t^\\ < c\\M l ^\\ 2 . The same 

estimate is obtained for \K 2 \ in a similar way. This shows (92) and completes the 
proof. ■ 

A. 2 Proof of Theorem 2.4 

For a fixed eigenvalue e ^ of L , define the subsets of N: 



AC (i) 

Mr 

Mi 



= {j\E z -E,=e}, 

= {i\Ei - Ej = e}, 

= Uj A/" r w = {jl-Ej - Ej = e for some «}, 

= Uf Af/ j) = {i\Ei - Ej = e for some j}. 
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We also let Pi denote the rank-one projector onto Ctpi, where we recall that {y?i} is 
the ortho normal basis diagonalizing H p . For any nonempty subset M C N, put 

is empty. 

jeAf 

Set E mn := E m — E n , and for e G a(L p )\{0}, m G Mi and n G A/" r , define: 

S m := Ma(P M ( m) GP K cGP M ( m) \P M ( m) ) > 0, (95) 
^ := Mafp^GPMcGP^ \P M ( n) ) >0. (96) 

Here, the superscript c denotes the complement. Notice that if e = 0, then M r c = Mf 
are empty, and 5 m ,5' n = 0. We define also 5 := inf meA ^{5 m } + inf neA f r {5' n }. From 
P(L P = e) = -. E . , =e y Pi <S> Pj, we obtain together with the definition of T(e) given 
in (16): 

m,n {i,j:E i:j =e} {k,l:E kl =e} 

The idea here is to get a lower bound on the sum over (m, n) G N x N by summing 
only over a convenient subset of N x N (notice that every term in the sum is positive). 
That subset is chosen such that the summands reduce to simpler expressions. 
Using the definition of m (see (17)), we obtain 

Pijin* P mn mP k i 
= Pij (Gig 1 — G r g 2 ) P mn (Gigi — G r g 2 ) Pm 
= PiGP m GPk <8> P n SjnSni\9i\ 2 — PiGP m <E> P n CGCPiSj n S m kgig2 

-P m GP k <g> PjCGCPJimSmg^ + P m <g> P 3 CGCP n CGCPA m 5 mk \g 2 \ 2 . 

Summing over i,j and k,l according to (97) yields 
5^ S Pijin* P mn mP k i 

{i,j:E lJ =e} {k,l:E kt =e} 

= (g^^GPm ®P n - g 2 P m <S> P N {m)CGCP^j ■ adjoint. 

For (m, n) e Mi x Mf, we have P N (n) = and P^w 7^ 0, and for (m, n) G Mf x A/" r , 
we have P^(n) 7^ and -P^-(m) = 0. As explained above, we now get a lower bound on 
the sum (97) by summing only over the disjoint union 

(m, n) G Mi X AC C U A^ c x M r . 
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An easy calculation shows that 



r p (e) > inf (f dS\g 2 (E t3 ,a)\ 2 ) ^ P m ®CP^ m) G P M * GP^ m) C 

+ inf (f dS\g 1 (E lJ ,a)\ 2 )Y, P ^ GP MfGP M < n) ®P n . 
E ar° \Js 2 J nlzK f 1 1 



neAf, 

Next, we investigate the integrals. From (10), we have 

/ dS\g 1:2 (E tJ ,a)\ 2 > \E l3 \ [ dS\g(\E i:j \,a)\ 2 , 
Js 2 Js 2 

uniformly in (5 > 1. With (95), (96) and remarking that a(CTC) = cr(T) for any 
selfadjoint T, this yields 

T p (e) > inf (\E tJ \ [ dS\g(E tJ ,a)\ 2 ) ( mf{5 m } + inf R}) P(L P = e), 

since EmeAf, P ™-® P m < T ) = UneMr p M W® P n = P(L p = e). This shows 1) of Theorem 
4.4. 

Now we look at the zero eigenvalue. A general normalized element of RanP(L p = 
0) is of the form = J2i c i¥i ® V 9 ^ with J2i \ c i\ 2 = 1> so 

S(E mn + u) (<pi <g> ifi, m*P mn mip j <g> ipj) . 



h3 



Using again the explicit form of m given in (17) and ((p m ,CGCip n ) = ((p m ,Gtp n ), we 
obtain 



r (°)0) = ( X ~~ 5e ^< °) y + M ) K^n^Vm)! 2 |C„5-1 - C m £f 2 | 



(9* 



We split the domain of integration R x S 2 into R + x S 2 U R_ x S 2 and using (10) 
and g 2 {u,a) = —gi(—u,a), arrive at 

J S(E mn + u)\cngi - c rn g 2 \ 2 

= J ^ 1 5 (E mn + u) a/ 1 + fic n g - y/JIc m g + 5(E mn - u) ^/Jlc n g - y/l + fic m g J . 

This together with (98) gives 



(0,r(O)0) = 2 Yl \(<Pn>G<Pr, 



|2 e 



{m,n:E mn <0} 



x |e _/3Em/2 c„ -e 
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- pEn/2 c m \ 2 j 5(E mn + co)\g\ 2 , (99) 



where we used S(E mn + u)fi = S(E mn + uj){e l3Emn — 1) 1 . Equation (99) shows that 
if we choose c n = Z p 1 / 2 e -P E n/2 ; then each term in the sum is zero. Recall now that 
the particle Gibbs state is given by (21), so (fijg, r(O)f^) = 0. Since T(0) > 0, this 
implies that Q 1 ^ is a zero eigenvector of T(0). 

Finally we show that there is a gap in the spectrum of T(0) at zero. Indeed, from 
(99), we get by the definition of g (see statement of Theorem 4.4): 



(0,T(O)0) > 2g Yl 

{m,n:E mn <0} 



-f3E m /2 _ -/3E n /2 I 



go 



m,n 



9o ^ ^ ^ -|- 6 ^ |c m | e ^ ^ (c n c m -|- C ra C m )) 

m,n 

= g (z p {(3) + Z p (/?) - 2\ J2^ Em/2 Cm\ 2 ) 

m 

= 2g Z p (f3)(l-\ 



where we used Yin l c «l 2 = 1- Therefore, we obtain on RanP^: T(0) > 2g Z p (P). 
This proves that if go > 0, then we have a gap at zero and zero is a simple eigenvalue. 



A. 3 Proof of Proposition 4.7 

We denote the spectrum of L p by cr(L p ) = {e,}, where we include multiplicities, i.e. 
for degenerate eigenvalues, we have tj = for different j ^ k. Let Pj denote the 
rank one projector onto span{(/?j}, where tpj G Tip <E> 7~C P is the unique eigenvector 
corresponding to ej. Let e be a fixed eigenvalue of L p . Setting rrij = Pj-rn, we have 



(100) 



Qi J m* ((L p -e + uf + e 2 ) 1 m Qitp 
= (^'Gi / m *j m j(( e j -e + u) 2 + t 2 Y 1 Q l iA . 

e,ea(L p ) \ J I 

First, we estimate the term in the sum coming from {j : Cj = e}: 

U,Qi j m* mj {u 2 + e 2 )- 1 Q l ^\ < Ju- 2 \\mjQM 2 . (101) 

{ej=e} ' {ej=e} 



Now 



W^QM 2 = \\P{L P = e)(G Wl - G r g 2 )Q^\\ 2 < 2||G|| 2 (M 2 + \g 2 \ 

{ej=e} 
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so (101) < 2\\Gnn 2 {\\gi/u\\l, + \\g2/u\\l 2 ) = 4\\G\\ 2 \\ gi /u\\ 2 L2 M\ 2 . From our 
assumptions on g (see (7)) and (10), it is clear that ||gi/-u|| L 2 = C < oo, uniformly in 
j3 > 1, and we conclude that 

(101)<C|H| 2 . (102) 
Next, we estimate the sum of the terms in (100) with ej ^ e and write it as 

/ du(( e i -e + u) 2 + e 2 )- 1 m i (M,V), (103) 

where we put rhj(u,ip) = f s2 dS\ \m,j(u, a)Qiip\ | 2 . V£ > 0, we have 
y / du ((ej — e + u) 2 + e 2 )~ 1 rhj(u, tp) 

rj/c -h\<>-(r-c,)\>Z} 



< C 2 J ||r^(^, q:)^!^!! 2 < 4e- 2 ||C|| 2 ||^ 1 /^||i 2 ||^M 2 < Ce- 2 ||^H 2 - (104) 
Next, with the changes of variables y — u — (e — e_j), we arrive at 
y / du ((ej — e + u) 2 + e 2 ) _1 m ,,•(«, tp) 

= (fdy (y 2 + e 2 )~ l \ ^m^e-e,,^) 

+ J dy ( y 2 + e 2 )- 1 ^ [fn j (y + e -e j ,4>)-m j (e-e j ,4>)]. (105) 



£ e^e 

The mean value theorem yields for the last sum: 



V d v\ye(-W J^^jiv + e- e j^)- ( 106 ) 



Now 



= 2 V / dS Re{Pj(d u m)(y + e - ej,a)Q 1 i/j,P j m(y + e ~ ej,a)Q!tp) . 
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Using the Schwarz inequality for sums, we bound the modulus of the r.h.s. from 

above by 



2 I dS jY J \\PAdum){y + e-e 3 ,a)QM\lY,\\ P My + ^-^^)QiW- 

(107) 

Now m(y + e — ej, a) = Gigi(y + e — ej, a) — G r g 2 (y + e — ej, a), so 

\\p im (y + e - e ^ a )QM 2 ( 108 ) 

< 2\ 9l (y + e - e v a^WP^Q^W 2 + 2\g 2 (y + e - e v a^WPjG^W 2 . 

We have to evaluate this at y = y G (— £,£)■ Clearly, |e — ej + y\ > |e — e 3 -| — |y| > 
rf — £ > if we choose £ < <i /2, where 

<i := inf | — e^- 1 > 0. 

The r.h.s. of (108) can thus be estimated from above by 

2 sup \g 1 (u,a)\ 2 \\P j G l Q 1 rl>\\ 2 + 2 sup \g 2 (u, a^W^GMW 2 , 

\u\>d /2 \u\>d /2 

hence we arrive at 

|(106)| < 32|y| ||G|| 2 |H| 2 / dS ( sup \d u9l \ + sup \ gi \) . (109) 

JS 2 \\u\>d /2 \u\>d a /2 J 

Using the conditions (7) with p > 0, one shows that the suprema are bounded, 
uniformly in (5 > 1, and so is \gi\, thus (109) gives 

|(106)| <C|y||M| 2 . (HO) 

Remark that the constant here depends on d , C ~ d p 1 ^ 2 . This argument is valid 
for any p. Going back to the second term on the r.h.s. of (105), we have shown: 

/ ,/l f 2 [™j(y + e ~ e J> ^) ~ ™j( e - e i> ^)] 

< C\W f J^ 2 dy<C%W. (Ill) 
Now we consider the first term on the r.h.s. of (105). We see that, as e/£ — > 0, 
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for any < rj < 1. This simply follows from the fact that for any such r], we have 
lim^oo x^Arctan^) — 7r/2) = 0. Also, 

^ mj(e - ej,ip) = J {ip, Qim*5(u - e + L P )P(L P ^ e)mQiip) . 
We conclude that (103) is equal to 

Choose e.g. £ = e 1//4 and rj close to 1, then the we arrive at 

(103) = 7 Qim*<S(« - e + L p )P(L p ^ e)mQ^) - 0(^/4) j |^j p | 

This together with (102) yields 

Qi y m*((L p -e + w) 2 + e 2 ) _1 m(5i 

> Qi^|ym*P(L p ^e)5(L p -e + M )m-0(e 1/4 )|g 1 . ■ 
A. 4 Operator calculus 

We outline an operator calculus for functions of selfadjoint operators, used extensively 
in this work. For a detailed exposition and more references, we refer to [HS3]. 

Let / G Co(R), k > 2, and define the compactly supported complex measure 
= —^(d x + id y )f(z)dxdy, where z = x + iy and / is an almost analytic 
complex extension of / in the sense that (d x + id y ) f(z) — 0, z G R. Then, for a 
selfadjoint operator A, one shows that 



f{A) = J df{z){A-z)-\ 



where the integral is absolutely convergent. Given /, one can construct explicitely an 
almost analytic extension / supported in a complex neighbourhood of the support of 
/. One shows that for p < k — 2, 

f \df{z)\ \lmz\-v- 1 < CJ2 \\f { %-P-i, (H2) 

where ||/|| n = / dx(x) n \f(x)\, and (x) = (1 + x 2 ) 1 ^ 2 . Furthermore, the derivatives of 
f(A) are given by 

f<?\A) =p\ f df(A)(A - z)-v-\ (113) 



We finish this outline by mentioning that these results extend by a limiting argument 
to functions / that do not have compact support, as long as the norms in the r.h.s. 
of (112) are finite. 
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